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ABSTRACT 

We calculated the influence of the limb-darkened finite disk correction factor in 
the theory of radiation-driven winds from massive stars. We solved the 1-D m-CAK 
hydrodynamical equation of rotating radiation-driven winds for all three known so- 
lutions, i.e., fast, f2-slow and 5-slow. We found that for the fast solution, the mass 
loss rate is increased by a factor ~ 10%, while the terminal velocity is reduced about 
10%, when compared with the solution using a finite disk correction factor from a uni- 
formly bright star. For the other two slow solutions the changes are almost negligible. 
Although, we found that the limb darkening has no effects on the wind momentum 
luminosity relationship, it would affect the calculation of synthetic line profiles and 
the derivation of accurate wind parameters. 

Subject headings: hydrodynamics — methods: analytical — stars: early-type — stars: 
mass-loss — stars: rotation — stars: winds, outflows 



-3- 



1. Introduction 



The CAK theory (ICastor et al] 1 19751) describes the mass loss due to radiation force in 
massive stars. This theory is based on a simple parameterization of the line force (a and k) 
which represents the contribu t ion of the spectral lines to the radiative acceleration by a power law 



distribution function. 



Abbott! (|1982|) improved this theory calculating the line force considering 



the contribution of the strengths of the hundreds of thousands of lines. He also included a third 
parameter (6) that takes into account the change in ionization throughout the wind. Despite this 
immense effort to give a more realistic representation of the l ine force, evident discr epancies 



still remained. Furthe r improvements to this theory done by 



Pauldrach et al. 



Friend & Ab bott (1986) and 



(1986) (hereafter m-CAK model) relaxed the point star approximation with the 
introduction of the finite disk correction factor, assuming a uniform bright spherical source of 
radiation. From then on, this model has succeeded in describing both, wind terminal velocities 
(Voo) and mass-loss rates (M) from very massive stars. As a result of the radiation force, the 
properties of the stellar winds must somehow refle ct the luminosities of the star s. This relationship 



can be obtained from the line driven wind theory (IKudritzki et al. 



1995, 



1999) and, nowadays, 



it is known as the Wind Momentum-Luminosity Relationship (WM -L). It predicts a strong 



dependence of wind momentum rate on the stellar luminosity with a (|Puls et al 



1996) 



The m-CAK hydrodynamical solution (hereafter the fast solution) is characterized by an 
exponential growth at the base of the wind that matches very quickly a /J-law profile when the 
velocity reaches some few kilometers per second, with a ft index in the range 0.8 to 1.0. 



However, in the last decade, 



Cure et al. 



(12011) found two new physical 



iCurel (120041) and 

solutions from the 1-D non-linear m-CAK hydrodynamics equation that describe the wind velocity 
profile and mass loss rates from rapidly rotating stars (the fi-slow solution) and from slowly 
rotating A- and late B-type supergiants (the 5-slow solution). The Q-slow solution only exists 
when the star's rotational speed is larger than ~ 3/4 of the breakup speed. This Q-slow solution 
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posses a larger mass loss rate (the higher the rotational speed, the higher the mass loss rate) and 
reaches a terminal velocity which is about 1/3 of the fast solution's terminal speed. On the other 
hand, the 5-slow solution is found when the line-force parameter 6 is slightly larger than ~ 0.25. 



Puis et al. 



iigh values of 8 are expected in hydrogen rich environments; for a pure hydrogen gas 
(|2000|) demonstrated that 6 is 1/3. This last solution, where the Abbott 6 factor represents changes 
in the ionization of the wind with distance, reaches a slow terminal velocity, similar to the Q-slow 
solution, but with a much lower mass loss rate. 

In the m-CAK model, the calculation of the radiation force is often carried out assuming a 
uniform bright finite-sized spheri cal star. The rapid rotation, however, c hanges the shape of the star 



to an oblat e configuration (ICranmer & Owocki 



darkening (|von Zeipel 



1995 



Pelupessy et al. 



2000) and induces gravity 



19241) as function of (co)-latitude. In both cases, a rotating and non-rotating 



star, the decrease of the temperature outwards the photosphere produces a limb darkening effect 
which also modifies the finite disk correction factor. The theoretical formalism for computing 
the self-consistent radiation force for non-spherical rotating stars, in cluding the effects of stellar 



1995) 



oblateness, limb darkening and gravity darkening, was developed by ICranmer & Owockil i 
However, to disentangle the effects of each one of these competing processes upon the wind 
structure, these authors present a semi-quantitative analysis and estimated that the limb darkening 
effect could increase the mass loss rate (M) in an amount of ~ 1 1% to ~ 13% over the uniformly 
bright models. However, th at larger mass loss would imply a reduction in the wind terminal speed. 
Owocki & ud-Doulal (|2004|) carried out (for the fast solution) a perturbation analysis of the effects 
of the gas pressure on the mass loss rate and wind terminal velocity in terms of the ratio of sound 
speed to escape speed (a/v esc ). They showed that for finite-disk-corrected spherical wind, typical 
increases in mass-loss rate are 10%-20%, with comparable relative decreases in the wind terminal 
speed. 



Then, considering that the radiative flux does not change significantly when the limb 
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darkening is taken into account, an enhancement of ~ 10% in the mass loss rate might lead 
not only to a lower terminal speed (v M ) but also to a change in the theoretical WM-L. An 
accurate determination of the WM-L relationship for A and B supergiants (Asgs and Bsgs) 
is important because it wou ld allow the use of these stars as extragalactic distance indicators 



( Bre solin & K udritzki 



2004) 



In this work, we present an analytical expression for the limb darkening finite disk correction 
factor and solve the 1-D hydrodynamical equation for all three known solutions for radiation 
driven winds; i.e., fast, tl-slow and <5-slow solutions. These results are compared with the wind 
solutions computed with the finite disk correction factor assuming a uniform bright star, finding 
that the effects of the limb darkening are only important for fast solution. 

In §2 we briefly describe the 1-D momentum equation of the wind, in §3 we present an 
analytical expression for the limb-darkened finite disk correction factor and in §4 we solve 
numerically the hydrodynamics equations for model parameters corresponding to the fast, Q-slow, 
5-slow and Q£-slow solutions. Finally, in §5, we discuss the results, conclusions and future work. 



2. The m-CAK hydrodynamic model 



The m-CAK model for radiation driven winds considers one dimensional component 
isothermal fluid in a stationary regime with sp herical symmetry. Neglecting the effects of 



viscosity, heat conduction and magnetic fields (ICastor et al 
conservation and radial momentum read: 



1975|), the equations of mass 



An r p v = M, 



(1) 



and 
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dv ldp GM(\-T) vj(r) , 

v — = — + + g ,me (p, dv/dr, n E ). 

dr p dr r L r 



(2) 



Here v is the fluid velocit y and dv/dr its gradient. All other variables have their standard 



meaning (see 



Cure (2004) for a detailed derivation and definitions of variables, constants anc 



functions). We adopted the sta n dard parametrization fo r the line force term, given by 



:1982), 



Friend & Abbott 



Abbott 



1986), 



Pauldrach et al. 



(119861) : 



f" e = £f»(r,v,dv/dr) [r 2 v^- 



n E 
W(r) 



(3) 



where the coefficient C depends on M, W(r) is the dilution factor and / D is the finite disk 
correction factor. 



Introducing the following change of variables u = —R*/r, w = v/a and W = dw/du, where 
a is the isothermal sound speed and a mt = v mt /a, where v rot is the equatorial rotation speed at the 
stellar surface, the momentum equation becomes: 

(1 \ d 2 / d \ a 

l-- I \w-£+A + -+al t u-C' f Bg (u)(wr d iw — \ =0. (4) 

The standard method for solving this non-linear differential equation © together with 
the constant C'(M) (eigenvalue of this problem) is imposing that the solution passes through a 
singular (or critical) point. 

Critical points are defined at the roots of the singularity condition, namely: 

d 

— F(u,w,w') = 0. (5) 
ow' 

At this specific point and in order to find a physical wind solution, a regularity condition must be 
also imposed, i.e., 

d OF dF , 

— F(u,w,w ) = — + — w = 0. (6) 
du ou ow 
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In order to solve this equation we need to know the behaviour of the finite disk correction 
factor / D . To disentangle limb darkening from rotational effects (gravity-darkening and 
oblateness) we will analize them independently. A discussion on the effects of th e oblate finite 



ArayaetaL 



disk c orrection factor on the velocity profile and mass loss rate was presented by 
(|201l|) . Therefore, in this work we mainly discuss the importance of the limb darkening on 
radiation driven winds. 



3. Limb-darkened finite disk correction factor 



Cranmer & Owockil (|1995|) derived an integral expression for the limb-darkened finite disk 
correction factor (/ld)> based in a simple linear gray atmosphere, namely: 

J2 



fu3 (r, v, dv/dr) = — 



(1 + aix' 2 f x 



1 + 2 



3 /yU' 2 -^ 



where cr = (d\nv/d\nr) - 1 and = Jl - R 2 /r 2 . 

The integration of eq. [7J gives the following analytical expression : 

1 / o-+ 1 



fx'dfi', 



Ad (r, v, dv/dr) = 



2cr(a+ 1) 



x 



A + cr+ 1 



,35 cr 

£r(ar+ 1) 2 F l \-,-a,-,- 

2 2 r l (cr + 1) - cr 



+r 2 (cr+ 1) 



where 2 F\ is the Gauss Hypergeometric function. 



cr+ 1 



-5 + cr+ 1 



- 1 +o- 



(7) 



(8) 



Figured] compares the run of both uniformly bright and limb-darkened finite disk correction 
factors as function u, using two different /Maw velocity profiles (J3 = 0.8 and 2.5) and a typical 
value of the a line-force parameter equals 0.6. This figure clearly shows that at the base of the 
wind the factor / LD (shown in gray-dashed line) is about ~ 10% larger than the one obtained 
for a uniform bright stellar disk, / D (continuous line), increasing, the value of the mass loss 
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rate. Instead, at larger distances from the stellar surface, both correction factors have the same 
behaviour as a function of u. 




\l=-RJl a=-RJi 



Fig. 1. — Uniformly bright fo (continuous line) and limb-darkened /ld (gray-dashed line) finite disk 
correction factors. Calculations where performed using a jS-law velocity field with /? = 0.8 (left panel) and 
P = 2.5 (right panel). In both cases the value for a was 0.6 



Mathematically, / LD can be considered as a small perturbation of / D , as it is shown in Figure 
[H e ven for different ff-law i ndex. Thus, based on the standard theory of dynamical system (see, 



e.g., 



J. Palis & de Melo 



1982|) . we expect no large differences when considering velocity profiles 
from the equation of motion (eq. SJ), for the cases where uniformly bright or limb-darkened 
finite disk correction factors are used. This is a consequence of the theorem of the continuous 
dependence of the solu tions of the ordinary differential equations (ODE) on their parameters 



(Hirsch & Smale 



1974) 



Therefore, the scope of this paper is limited to the study of the numerical 1-D stationary 
solutions of eq. |4] leaving a theoretical topological analysis (and also a time-dependent one) for a 
future work. 
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4. Results 

We are now in conditions to solve the non-linear differential equation (eq. ID) considering 
the factor / LD given by eq. [8j The calculation of all partial derivatives of / LD (w, w, w') are given 
in Appendix [Aj These derivatives are needed in order to evaluate the singularity and regularity 
conditions (Eqs. [5]and[6l respectively). Depending on the selected parameter-space, each type of 
solution explains the wind of a different kind of massive object, i.e., fast solution describes the 
wind of hot stars, the tl-slow solution explains the wind of rapid rotators such as Be stars, and the 
5-slow solution characterizes the wind of A-type supergiants. In the following subsections we will 
adopt a prototype star for each one of these three known physical solutions, in order to analyse the 
effects of the limb darkening on the m-CAK hydrodynamical model. 



4.1. Fast Solution 



For the standard fast solution we selected, as in 



Curd (12004 . a typical 05 V star with the 



following stellar and line-force para meters: r sff = 45 000 K, 



k = 0.124, a = 0.64 and 6 = 0.07 (Lamers & Cassinelli 



to solve the momentum equation is described in 



og g = 4.0, R/R Q = 12, v rot = 0, 



1999). The numerical code we used 



Curd (|2004l) . Figure [2] (left panel) shows the 



velocity profile for the standard case, where a uniform bright star disk (continuous line) and the 
limb-darkened one (gray-dashed line) are used. Figure [2] (right panel) displays the difference in the 
velocity, Av = v un - v LD (where v un is the velocity profile when the uniform finite disk correction 
factor is taken into account, while vld is the wind solution obtained using /ld)- The vld(w) profile 
is always smaller than the v un (u) one, with a monotonically increasing difference. The effect of 
the limb-darkened finite disk correction factor changes the behaviour of the velocity field in the 
most external layers, it reaches a smaller terminal velocity by about 10% of the v M value of the 
standard m-CAK case. There is no significant change of the velocity field at the base of the wind. 
Therefore, the location of the singular point is almost the same in both cases. Our calculations 
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confirmed the predictions of ICranmer & Owockil (|1995|) : lOwocki & ud-Doulal (120041) based on 
the behaviour of the / LD at the base of the wind, i.e., that the mass loss rate is increased a factor of 



about ~ 10%. Concerning to the WM-L relationship, the value of D mom = {Mv^ y[RJRo) shows 
almost no change due to a compensation of the increase in the mass loss rate and the decrease of 
the terminal velocity, as shown in Tabled! Although the value of D mom seems to remain unaltered, 
we would expect minor differences in the synthetic spectra when they are computed with the two 
different velocity profiles. 





Fig. 2. — Left panel: velocity profile as function of the inverse radial coordinate u. The standard 
m-CAK model is shown in continuous line and the solution with the limb-darkened finite disk cor- 
rection factor is in gray dashed line. The effect of the /ld in the velocity profile is very significant 
reducing the terminal velocity approximately in 10% with respect to the standard m-CAK model. 
Right panel: Velocity difference between the standard solutions with a uniformly bright, / D , and 
the limb-darkened, / LD , correction factors. 



These results show that the correction to the line radiation force due to the limb darkening 
effect leads to lower mass loss rates and higher wind terminal velocities, both in approximately 
10%, when compare with the contribution of a uniformly bright star disk radiation source. 
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Table 1 : Wind parameters for the fast solution with uniformly bright (/d) and limb-darkened (/ld) finite 
disk correction factors 



/d /ld 



M(10- 6 M a yr- 1 ) 


2.206 


2.449 


Voo (km s ) 


3384 


3071 


'"singular 


1.027 


1.031 


Eigenvalue (C) 


40.89 


38.53 


log D mom (cgs) 


29.21 


29.22 



4.2. Q-Slow Solution 



The Q-slow solution is present when the star is rotating at velocities near the breakup 
rotational speed. Therefore, to study the effects of the limb darkening in the radiation force we 
select, the case of a typical Bl V star with high rotational speed (Q = v rot /vb r eakup = 0.9) and the 
following stellar parameters: T eS = 25 000 K, log g = 4.03, R/R p = 5.3. The c orresponding line 



force parameters: k = 0.3, a = 0.5 and 6 = 0.07 were taken from 



Abbott! (119821) 





Fig. 3. — Left panel: Same as Figured] v(«) versus u. In this case the effect of the / LD in the 
velocity profile is minimal. Right panel: Velocity difference. 



The resulting velocity profiles with uniform and limb-darkened correction factors and the 
corresponding differences in the velocities are shown in Figure [31 These plots show clearly that 
the effect of /ld in the velocity profile is minima. The influence of the / LD on the mass loss rate 
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Table 2: Wind parameters for the Q-slow solution with uniformly bright (/*d) and limb-darkened (/ld) 
finite disk correction factors 





/d 


/ld 


M(1(T 6 M a yr- 1 ) 


4.22 1(T 3 


4.22 10" 


Voa (km s~ 1 ) 


446.8 


446.5 


'"singular (^*) 


26.14 


26.14 


Eigenvalue (C) 


78.31 


78.27 


log D mom (cgs) 


25.44 


25.44 



and other wind quantities are shown in Table El together with the comparison of the velocity 
profile using the uniform correction factor. All the changes in these quantities are minimal or 
even negligible. There is an important dominance of the centrifugal force term in the momentum 
equation ©. 



4.3. 5-Slow Solution 

For the calculation of the /ld correction factor in the parameter-space of the 5-slow solution, 
we select an A-type supergiant star with the following fundamental parameters: = 10000 K, 
log g = 2. 0, R/Rn = 60, = 0, and line-force parameters: k = 0.37, a = 0.49 and 8 =0.3 (model 



W03 fromlCure et al. 



201 lh 



Similar to the Q-slow wind solution, the effect of the limb darkening is negligible in both, the 
velocity profile and mass loss rate (see Figured and Table [3]). 

Concerning the influence of the limb darkening on the WM-L relationship, there is no 
substantial effect. 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 



Fig. 4. — Left panel: Same as Figure [Q v(w) versus u. In this case the effect of the / LD in the 
velocity profile is minimal. Right panel: Velocity difference. 

Table 3: Wind parameters for the 5-slow solution with uniformly bright (/d) and limb-darkened (/ld) finite 
disk correction factors 





/d 


/ld 


M(10~ 6 M a yr- 1 ) 


7.22 ur 4 


7.36 10- 


Voa (km s~ 1 ) 


203 


200 


'"singular (^*) 


11.06 


11.06 


Eigenvalue (C) 


63.78 


63.54 


log D mom (cgs) 


24.85 


24.86 



4.4. £2£-Slow Solution 

Here we investigate the particular case when Q and 5 take higher values. We selected the 
same test star as in §!4.2l but with a different value of the 5 parameter (5 = 0.25). The computed 
hydrodynamic solutions for uniformly bright and limb-darkened correction factors are almost the 
same, as it shown in Figure \5\ and Table 0] 

When we compare the velocity profiles between the Q-slow solution computed with 6 = 0.07 
(see figure|4]left panel) and 8 = 0.25 (see figure [5] left panel), we find that both profiles have the 
same behaviour as a function of r. Therefore, the centrifugal term due to the high dominates over 
the 5-factor in g Une . Nevertheless, the influence of the ^-factor is not negligible, it reduces the mass 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 



Fig. 5. — Left panel: Same as Figured] v(w) versus u. In this case the effect of the / LD in the 
velocity profile is minimal. Right panel: Velocity difference. 

Table 4: Wind parameters for the £25-slow solution with uniformly bright (/d) and limb-darkened (/ld) 
finite disk correction factors 





/d 


/ld 


M(10~ 6 M yr- [ ) 


8.63 icr 4 


8.64 1Q- 


Voo (ims -1 ) 


367.8 


367.5 


r singular (^*) 


38.17 


38.17 


Eigenvalue (C) 


113.9 


113.9 


log D mom (cgs) 


24.66 


24.66 



loss rate in ~ 80% and the terminal velocity in ~ 20%. 

5. Discussion and Conclusions 

In this work we improved the description of the radiation force taking into account the 
correction factor due to a limb-darkened disk. In particular, we derived an analytical formula to 
compute this contribution. Then, we solved the 1-D non-linear momentum equation for radiation 
driven winds and analized the influence of /ld for all the three known solutions, namely: the fast, 
the Q-slow and the 5-slow solutions, as well as, the case of a high Q. and 8 parameter, the Q. 5-slow 
solution. 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

u=-RJl u=-RJr 

Fig. 6. — Left panel: Normalized velocity gradient dw/du versus u from the solutions of the 
equation of motion using the uniform finite disk f D . In continuous-line is plotted the gradient of 
the fast solution; dashed-line correspond to the £-slow solution; dotted-line to the tl-slow solution 
and dashed-dotted-line to the Q^-slow solution. Right panel: id, for the cases where the f LD us 
used in the equation of motion (eq. H]) 

We selected the appropriate stellar parameters of massive stars that are representative of 
each possible hydrodynamical solution and evaluated the velocity profile as function of the radial 
coordinate. 

We found a significant impact of /ld in the radiation driven-wind of massive stars that 
are described by the fast solution. Due to the effect of a limb-darkened disk the mass loss 
rate increased in an amount of ~ 10% while the terminal velocity is reduced about the same 
factor. Therefore, the limb darkening effect should be considered always in the calculation of the 
hydrodynamics fast solution. 

On the other hand, the influence of / LD on the fi-slow and 5-slow solutions is minimal. 
The maximum difference obtained in the velocity profile computed with uniformly bright and 
limb-darkened disk radiation sources is less than 3 km s~ l for the Q-slow solution, 7 kms~ l for the 
8- slow solution and 1.5 kins' 1 for the case when both parameter 6 and £1 are high (the Q£-slow 
solution). Therefore, the limb darkening effect is negligible when computing the wind parameters. 
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However, rotational effects like the star's ob 
wind in the polar direction (see lAraya et al 



ateness should be considered, since it modifies the 



201 if) being much faster than the spherical one. 



Moreover, the slow solutions predict even slower and denser flows than the spherical ones. 

The influence of / LD on radiation driven winds can be interpreted in terms of the resulting 
velocity profile. The mayor differences between the uniformly bright and limb-darkened finite 
disk correction factors are in the region just above the stellar photosphere, as Figured] shows. In 
this region the velocity from all the models described in section |4] are small, however the value 
of the velocity gradient from the fast- solution is 5 to 10 times larger than the values from any 
slow-solution. Figure [6] shows the normalized velocity gradient dw/du as function of u for the 
different types of solutions. Is this dependence on the velocity gradient, specifically in the finite 
disk correction factor, that makes a significant difference in the terminal velocity and the mass 
loss rate only for the fast solution and not for the slow ones. 

Concerning the WM-L relationship, the limb-darkened correction factor has no effects. 
The increase produced by the fast solution on M is compensated by a similar decrease of v^,. 
Considering the importance of having a theoretical WM-L relationship for B- and A- type 
supergiants the effect of the star's oblateness and gravity-darkening should be explored, together 
with the calculation of the synthetic line spectrum in order to derive accurate wind parameters. 
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A. Partial derivatives of f LD 

In order to find the location of the singular point, we need to evaluate the singularity 
condition given by eq. [5] and, then, impose the regularity condition given by eq. [6j To perfom this 
calculation we need to know all the partial derivatives of / LD (w, w, w'); i.e., df LD /du, df LD /dw 
and dfuoldW. 

Defining the following auxiliary variables: 

Z = w/w' (Al) 
A = u(u + Z) (A2) 

Thus, in terms of A, the finite disk correction factor for a uniformly bright spherical star / D reads, 

h&^Trr^il 1 -^-**™]' (A3) 

(1 + a) A 1 J 
while the limb-darkened finite disk correction factor / LD is: 

(1 - A) a [d - AT a + A(a + 1) 2 Fi (§, -a, |,^) + ^-ll 

/ld (A) = * ... ^ n V2 * (A4) 

2/l(a + 1) 

Defining now e (/I) = 5/ L d (A)/dA, we obtain: 

CW " 4( tt+ < l)a-l)^ X ' (2 - (3tt + 5) - i)( '-- <r + 



'3 5 /i 

+2 (A - 1) (a^ + 1) + (a + 1) A (2aA + 3) 2 Fj I -, -a, -, 



(A5) 



X 2A-1 

Therefore, all the partial derivatives can be calculated using the chain rule, getting, 

,,v 1 dfhD W dfui w' 2 dfu3 fKC . 

2m + w/w aw m aw uw ow 
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